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1. Introduction 


ABSTRACT 


Electrical polarization phenomena in GaN are important as they have significant impact on the operation 


of modern day energy efficient lighting and are fundamental to GaN-based high power and high fre- 
quency electronics. Controlling polarization is beneficial for the optimization of these applications. GaN 
is also piezoelectric, and therefore mechanical stress and strain are possible handles to control its po- 
larization. Nonetheless, polar semiconductors in general, and GaN in particular, are weak piezoelectric 
materials when compared to ceramics, and are therefore not considered for characteristic electrome- 
chanical applications such as sensing, actuation and mechanical energy harvesting. Here, we examine 
the effect of nanoscale porosity on the piezoelectricity of initially conductive GaN. We find that for 40% 
porosity, the previously conductive GaN layer becomes depleted, and exhibits enhanced piezoelectricity 
as measured using piezoresponse force microscopy, as well as by using a mechanical energy harvest- 
ing setup. The effective piezoelectric charge coefficient of the porous GaN, d33 ef, is found to be about 8 
pm/V which is 2-3 times larger than bulk GaN. A macroscale device comprising a porous GaN layer de- 
livered 100 nW/cm? across a resistive load under a 150 kPa mechanical excitation. We performed finite 
element simulations to analyze the evolution of the piezoelectric properties with porosity. The simula- 
tions suggest that increased mechanical compliance due to porosity gives rise to the observed enhanced 
piezoelectricity in GaN. Furthermore, the simulations show that for stress-based excitations, the porous 
GaN electromechanical figure of merit is increased by an order of magnitude and becomes comparable to 
that of barium titanate piezoceramics. In addition, considering the central role played by GaN in modern 
electronics and optoelectronics, our study validates a very promising research direction when consider- 
ing stress-based electromechanical applications which combine GaN’s semiconducting and piezoelectric 
properties. 


© 2020 Elsevier Ltd. All rights reserved. 


electric, and therefore mechanical strain or stress correspond to 
development of electrical polarization or field [4,5]. Interestingly, 


GaN and related III-N technology is becoming increasingly im- 
portant in electronic and optoelectronic applications, with emerg- 
ing high-frequency and high-power applications in addition to its 
established role in energy-efficient lighting [1,2]. Electrical polar- 
ization plays an important role in GaN technology, due to inherent 
properties of III-N crystal structure, and effects brought about by 
heterostructure growth |3]. In particular, III-N materials are piezo- 
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polarization fields in GaN based light emitting diodes (LEDs) are 
considered detrimental to LED performance [1,3], whereas they 
are a fundamentally enabling effect in GaN high electron-mobility 
transistor (HEMT) technology. The latter is implemented in GaN 
without the need for doping by making use of an interfacial po- 
larization discontinuity [3]. This indicates the importance and ben- 
efits of controlling this phenomenon. 

Furthermore, semiconductors which exhibit piezoelectricity 
make excellent smart/multifunctional materials, as they combine 
physical characteristics of distinct disciplines, and open up non- 
trivial routes to control these properties through their coupling. 
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For example, piezoelectric semiconductor nanowires (NWs), and 
other geometries, have been extensively studied in recent years 
as active components in applications where mechanical stress or 
strain controls the electronic/optoelectronic characteristics of semi- 
conductor NW devices (such as their current-voltage characteris- 
tics, or light emission intensity). These applications are referred to 
as piezotronic/photo-piezotronic |6-11]. III-Ns and non-nitride III- 
Vs are attracting increased attention in piezotronics [5,12] due to 
their mature processing, possible silicon integration and span of 
desired electronic and optical properties [13-15]. 

Porosification has been used as a processing tool to control the 
optical properties of GaN [16-19] as well as its mechanical prop- 
erties [20]. Furthermore, several reports demonstrated increased 
electromechanical coupling following GaN porosification [21-23]. 
In this report we examine the effect of nanoscale porosification 
on the piezoelectric properties of GaN, both at the nanoscale and 
macroscale, experimentally and through computational modelling. 
We discuss the results specifically in the context of mechanical to 
electrical energy transfer efficiency and compare the performance 
of porous GaN to conventional ferroelectric piezoceramics. Our re- 
sults suggest that porosification is an effective method to increase 
GaN efficiency as a piezoelectric material and that the performance 
could be expected to become comparable with ferroelectric piezo- 
ceramics, depending on the application. Considering the significant 
role of GaN in modern electronic technology, this finding has con- 
siderable implications for integrated GaN-based electromechanical 
technologies. 


2. Theoretical background 
2.1. Piezoelectricity and electromechanical coupling 


Piezoelectricity is the linear coupling of electrical polariza- 
tion/field (P/E) with mechanical stress/strain (T/S) fields [4]. Piezo- 
electricity is manifested in non-centrosymmetric crystals, with 21 
out of 32 crystal structures exhibiting piezoelectricity. The lin- 
ear relations between mechanical stress/strain (T/S) and electrical 
field/displacement (E/D), are given by 


T-YS (2) 


where gg is the vacuum permittivity and er the relative permittivity 
and Y Young’s modulus. The constitutive piezoelectric relations are 
given by [24] 
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where the ij indices identify the piezoelectric coefficient in the 
3 x 6 matrix representing the piezoelectric tensor (Voigt notation). 
Subscripts i {1-3} and j {1-6} represent electrical and mechani- 
cal fields correspondingly. The right-hand (left-hand) side equation 
refers to the converse (direct) piezoelectric effect, and superscripts 
refer to invariable fields. For example, the direct piezoelectric volt- 
age coefficient, gj, describes the change of the i component of the 
electric field due to a change in the j component of the mechanical 
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stress field, under constant or zero electrical displacement field - 
corresponding to open-circuit conditions, i.e., the piezoelectric as 
a voltage source. Similarly, the converse piezoelectric strain coeffi- 
cient (also called the piezoelectric charge coefficient), dij, describes 
the j component of the mechanical strain developed in response to 
the i component of an applied electric field, under constant or zero 
mechanical stress. 

Atomic force microscopy (AFM) can detect movement on pi- 
cometer order. The experimental conditions mentioned above cor- 
respond to piezoresponse force microscopy (PFM) operation, where 
the AFM tip is in contact mode, and follows the surface deforma- 
tion induced by the applied voltage through the tip. Therefore, PFM 
allows a direct measurement of the out-of-plane d coefficient of 
the converse piezoelectric effect. 

Piezoelectricity in GaN 

GaN usually crystallizes in wurtzite structure [3], which deter- 
mines the form of its piezoelectric matrix. The characteristic ma- 
trix for wurtzite structures takes the form 

0 0 0 0 —d;5 O 
d= 0 0 0 -di5 0 0 (7) 
—d3; —d3,; d33 O 0 0 


with experimental values for GaN that range |d3,| =1.5—1.9 
pm/V, d33 = 2.5 — 3.8 pm/V, |d45| = 18-3 pm/V [25-27]. 

As mentioned above, the combination of semiconducting and 
piezoelectric properties yields interesting piezotronic and photo- 
piezotronic phenomena. To some extent, this combination is some- 
what contradictory, as good conductors make poor piezoelec- 
tric materials, due to internal screening. Indeed, previous reports 
studying porous GaN piezoelectricity link it with carrier depletion 
[21-23]. 


3. Results 
3.1. Piezoresponse force microscopy 


Porosified and as-grown (non-porous) GaN samples were exam- 
ined using PFM. The porous layer was about 260 nm thick. Fig. 1a 
shows a cross section scanning electron microscopy (SEM) image of 
the porous GaN layer. Fig. 1b shows a schematic of the layer struc- 
ture connected for PFM: the only relevant layers are the remnant 
n-doped GaN (i.e. the part not etched), and the porous layer. The 
n-doped layer is electrically connected to the AFM stage (where 
voltage is applied) and the tip (electrical ground) is scanned along 
the porous layer. The porosity of this sample was estimated as 40%, 
see Supporting Information Section S1. 

PFM was carried out using 2, 4, 6, 8 V AC drive voltages, and 
the deflections across a 250 x 500 nm? area were recorded 
(see further details in Methods section). When performing PFM 
measurements, there are two major parasitic effects to consider: 
i) background - a parasitic offset signal coming from the instru- 
mentation |28]; ii) electrostatic contribution to the signal (arising 
from contact potential differences) |29]. Under the assumption that 
the background is mostly dependent on the experimental config- 
uration, we measured the sample before etching and subtracted 
the obtained results from the measurements. The as-grown (un- 
etched), n-doped, GaN layer is not expected to exhibit piezoelec- 
tricity due to screening associated with its high conductivity (car- 
rier density of 5-10'® cm~), as piezoelectricity sharply decreases 
with conductivity [30,31]. See more details in Supporting Informa- 
tion S2. We therefore use the following expression for d33 


d33 = Cmeasured = background (8) 


In order to minimize the electrostatic contribution, we carried 
out Kelvin probe force microscopy (KPFM) measurements prior to 
the PFM measurements, and used the result as a DC bias [29]. An- 
other consideration is signal calibration: we used a polyvinylidene 
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Fig. 1. (a) Cross section SEM image of the etched sample. Similar samples were used in the experimental parts of the paper. The red mark indicates 35 nm; (b) Schematic 
of PFM operation above the layer structure. The “undoped” layer is in practice unintentionally doped (UID) to a level of about 1-1016 cm~. 
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Fig. 2. PFM deflection vs. AC voltage for the porous GaN sample (blue data points), 
linearly fitted to yield 8 pm/V coefficient (slope of the blue dashed line). The range 
of theoretical bulk values is shown for reference (area). 


fluoride-trifluoroethylene (PVDF-TrFE) film made in house [32], un- 
der the assumption that after poling its piezoresponse is 22 pm/V 
[33,34]. See Supporting Information S2 for details. 

Fig. 2 shows calibrated and corrected PFM deflection vs. applied 
voltage. Each data point represents the spatial mean and stan- 
dard deviation of a single PFM scan, obtained in a resolution of 
256 x 512 points (scans are shown in Supporting Information Fig- 
ure S6). The piezoelectric coefficient was obtained by the slope of 
linear fitting of the data points. This yielded a porous GaN coeffi- 
cient of d33 =8+1.6 pm/V. The relatively large error margin was 
due to taking into account a 20% uncertainty in the d33 coefficient 
value of PVDF-IrFE used here for PFM calibration |32]. In any case, 
the standard deviations of the PFM measurements are smaller than 
20%. The corresponding curve for bulk GaN (filled area calculated 
based on known coefficient range) is also shown for reference. Our 
results suggest that for ~40% porosity, the effective d33 is about 
three times that of bulk GaN. 


3.2. Mechanical energy harvesting 


We further examined the macroscale electromechanical perfor- 
mance of the etched GaN layer. Current-voltage (J-V) measure- 
ments of the sample did not show any measurable current (data 
not shown), and therefore piezotronic activity [10] was not rele- 
vant for this sample. Considering that the starting material was a 
5 - 10'® cm~? Si doped GaN layer, this result indicates the layer 
became depleted. 

Fig. 3 shows the short-circuit current measured while the layer 
was periodically tapped by a blunt indenter applying 1.5 N at 2 Hz. 
The indenter was roughly the size of the contact (2 x 5 mm7), 
leading to an assessment of the applied stress of 150 kPa. Notice- 
ably, the approach and retract (compress and release) peaks are 
not symmetrical. The experiment was repeated with a set of resis- 
tive loads spanning 1 kQ2 - 100 MQ. A peak-finding algorithm was 
used to locate the highest approach and retract peaks (as seen in 


Fig. 3b), and the RMS value was calculated accordingly. Although 
the mechanical excitation was applied at 2 Hz, only the actual du- 
ration of contact is relevant for power transfer, and so the signal 
was integrated around the two prominent peaks to find RMS cur- 
rent. 

The maximum power transfer from the 150 kPa stress is about 
9 nW, obtained for loads of 5 and 10 MQ (indicating the optimized 
value would be slightly higher). Considering the contact area, this 
corresponds to an area density of 100 nW/cm? and a volume den- 
sity of about 4 mW/cm?, accounting for about 260 nm of porous 
layer thickness. These results (for area density) are in good agree- 
ment with previous porous GaN energy harvesting experiments re- 
ported in recent years [22,23]. 


3.3. Finite element simulation - COMSOL 


In order to obtain a deeper understanding of the electrome- 
chanical effects at play in porous GaN, we carried out COMSOL 
simulations on GaN structures with varying porosity. As shown in 
Fig. 4, we simulated a hexagonal unit cell (100 nm long diago- 
nal dimension), with periodic boundary conditions - effectively a 
thin film. The pore was cylindrical with a varying diameter, en- 
abling the presentation of results as a function of porosity, which 
is calculated as the volumetric ratio of air to GaN. This design is 
inspired by the observed morphology (Fig. 2). The total thickness 
was 500 nm and the pore depth was 400 nm. The structure was 
capped by a 10 nm GaN layer to simplify the application of forces 
to the structure (constant stress). The pore diameter was varied be- 
tween 1 and 81 nm, corresponding to porosity of up to 42%. The 
default GaN parameters were used with the addition of the piezo- 
electric coefficients d3; = 1.5, d33 = 3.8, djs = 1.8 pm/V (with the 
appropriate sign applied). A relatively high d33 value was chosen 
to increase the simulated effect. Fig. 4b shows the potential distri- 
bution for the largest pore size. 

Fig. 5 shows the ‘raw’ output of the simulation: the top sur- 
face displacement and potential difference across the simulated 
element as a function of porosity - for a given load of 5 GPa. 
The results are straightforward - deflection and voltage both in- 
crease considerably with porosity (two, and three-fold at 42%, cor- 
respondingly). This would be expected, as a porous material is 
more compliant [35,36], and the remaining material produces the 
open-circuit voltage. This raises several interesting questions, for 
example what are the effective piezoelectric constants? and what 
is the electromechanical coupling efficiency as a function of poros- 
ity? 

To examine these questions, we calculated several quantities 
based on the simulation results. First, we consider the effective 
out-of-plane piezoelectric coefficient. The electrical boundary con- 
ditions for the model were open-circuit, therefore calculating d co- 
efficients is not straightforward and requires calculating the piezo- 
electric voltage coefficient, g, first (see Eq. (5)), since the d coef- 
ficient is defined under constant electric field. We have recently 
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Fig. 3. Generated current measured during repetitive loading of the porous GaN layer, by applying 1.5 N / 150 kPa: (a) measured short-circuit current. The interrupted 
intervals are due to data recording and storage, the excitation was continuous; (b) zoom into one of the peak regions showing compress and release peaks. The prominent 
peaks were identified and used to calculate RMS current as well as the power; (c) RMS current (blue dots, left axis) as a function of load resistance and the corresponding 
power (orange asterisk, right axis); (d) top view of the device, where the left hand side arrow (red) points to the sputtered top contact and force application location, and 
the right hand side arrow (yellow) points to the bottom contact - effectively touching the un-etched part of the sample. 
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Fig. 4. COMSOL simulations: (a) the simulated geometry, a periodic GaN structure, with an embedded pore of varying diameter. The arrow shows the direction of 5 GPa 
stress applied to the structure; (b) the piezoelectric potential developed following the applied pressure for the largest porosity value of ~42%. 


demonstrated this procedure experimentally [37]. The g coefficient 
is the ratio between the developed electric field and the applied 
stress. The relation between the two coefficients is given by d = €g, 
and this allows us to calculate d, as shown in Fig. 6. The piezo- 
electric coefficient follows the trend seen in the displacement and 
piezoelectric voltage - an increase with porosity. We point out to- 
wards two interesting observations - the highest value found is 
7.93 pm/V for 42% porosity, in very good agreement with the ex- 
perimental value found using PFM for about 40% porosity. Also, the 
zero porosity value, ~ 3.4 pm/V, is actually lower than the input d33 


(3.8 pm/V). As discussed below in detail, this is a demonstration of 
the bulk confinement effect, where the lack of it contributes to the 
increase in piezoelectric coefficients of nanomaterials [38]. 

Next, we examine the energy conversion efficiency. We consider 
two figures of merit for the electromechanical (EM) coupling: one 
is the electromechanical coupling efficiency coefficient, k?, defined 
as the ratio between stored electrical energy and input mechanical 
energy; the other is the stress-based energy harvesting figure of 
merit (EH FoM) which is given by the product d - g [39,40]. These 
quantities take the form 
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Fig. 5. Simulated effective d33 as a function of structure porosity. 
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Fig. 6. Simulated displacement (right hand side) and piezoelectric potential (left 
hand side) in the studied structure as a function of volume porosity for 5 GPa stress. 
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Fig. 7. shows the two figures of merit calculated based on the simulation, as a func- 
tion of porosity. For the efficiency calculation, we considered the input mechanical 
energy as the product of the applied force and the displacement, and the stored 
electrostatic energy was calculated directly by COMSOL; for FoM?. calculation, we 
approximated the effective permittivity as a linear function of porosity volume frac- 
tion (7), Ee¢r = NEcan€o, and used Equation 10. 





d?. 
2 1J 
T Eisi 
2 
FoM;, = dij '8&ij = zT (10) 


ii 
Noticeably, the difference between the two figures of merit 
is the lack of the compliance term (si) in the denominator in 


FoM}, which means that softer materials do not suffer an effi- 


ciency penalty, making compliant materials such as PVDF good 
sensors/energy harvesters [39,41,42] - when the excitation is based 
on stress and not strain. 
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The results indicate a significant distinction between the two 
figures of merit: while the electromechanical coupling efficiency 
remains nearly unchanged with porosity (slight increase), FoM}; in- 
creases nearly 10 fold. This suggests that porosification is a promis- 
ing strategy for increasing efficiency in stress based piezoelectric 
applications. 


4. Discussion 
4.1. Transition of GaN from conductive to an active piezoelectric 


The porosification process described here is based on conduc- 
tivity, and its efficiency increases with doping/conductivity. We 
consider two routes for the transition from a conducting starting 
material to an un-screened piezoelectric material: i) the remain- 
ing material has increased effective resistance due to geometry (di- 
mensions or connectivity) ii) the remaining nanostructure becomes 
depleted, due to increased surface area (and hence surface deple- 
tion) and possibly due to size effects (quantum confinement). We 
have observed that similar etching profiles obtained from samples 
of initially higher doping levels then reported here, result in struc- 
tures that are still conducting. This is an indication that the geom- 
etry and connectivity aspects alone, do not explain the significant 
decrease in conductivity and onset of piezoelectricity. We therefore 
focus on the second possibility, size induced depletion. 

The onset of carrier depletion in GaN due to porosification has 
been considered to be an important mechanism in previous work 
regarding porous GaN piezoelectricity [23]. The relevant feature 
size in our material, as seen in Fig. 1a, is about 35 nm, although 
some features are thinner than this. Previous studies have shown 
that this size falls in the range of critical thickness (in form of GaN 
nanowire diameter) resulting in a significant decrease in conduc- 
tivity due to surface depletion [43,44]. This is a general outcome of 
size and geometry in semiconductors and not unique to GaN [45]. 
Therefore, we agree with previous interpretations that size related 
carrier depletion is responsible for the onset of piezoelectricity in 
the initially doped GaN. In the following section, we discuss the 
emergence of piezoelectricity in porous GaN which is considerably 
larger than bulk GaN. 


4.2. Increased piezoresponse of porous GaN 


Lateral mechanical clamping reduces the ability of a piezo- 
electric thin film to react to an excitation. It has been demon- 
strated that nanostructures experience reduced clamping, and an 
increased piezoelectric response [38]. Lateral clamping is expressed 
by [46] 


2513 qMatrix (9) 


Eff _ 4Matrix 
d33 d i 31 
S11 + S12 


33 


where s; are compliance matrix elements, and superscripts indi- 
cate matrix elements and effective/measured values. The matrix 
elements we used in our COMSOL simulation are do = 3.8, 
aa x =-1.5 pm/V and sz, = 3.03, S42 = —0.659, s43 = —0.322 
GPa~!. Using Eq. (9), this yields dil! = 3.39 pm/V, which is in 
very good agreement to the simulation result (Fig. 6 at zero poros- 
ity), indicating that the simulation captures this effect successfully. 
Logic dictates that substrate clamping removal will allow the mea- 
sured/simulated coefficient to be only as high as the matrix coeffi- 
cient (3.8 pm/V in our case). Therefore, it does not explain the full 
extent of the results as they are manifested in PFM and simulation. 
The conclusion therefore is that the overall structure of the porous 
material, and the increased compliance in particular, are responsi- 
ble for the results. 
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4.3. Comparison to porous ceramics 


As mentioned above, porosity as a tool to enhance piezoelec- 
tricity has been studied mostly in the context of ceramics |40,47- 
49], polymers [50], and even as a route to induce piezoelectric- 
ity in originally centrosymmetric 2D materials [51,52] - in addi- 
tion to GaN [22,23]. When considering ceramics, and ferroelectrics 
in general, a poling process is necessary for dipole alignment 
and onset of measurable piezoelectricity. This process is depen- 
dent upon the ability to efficiently apply an electric field across 
the ferroelectric material. The ability to efficiently pole a porous 
structure is found to depend both on the total porosity, and the 
pore shape, where aligned/directional pores offer better efficiency 
in poling [47,49]. This then affects the piezoelectric properties of 
these structures, such that the increased compliance (and subse- 
quent increased piezoelectricity) competes with the reduced poling 
efficiency, brought by inefficient electric field distribution [40,48]. 

Conversely, GaN is not a ferroelectric, and once grown, it does 
not require poling. Therefore, the only effect at play is the in- 
crease in mechanical compliance. As seen above, the electrome- 
chanical coupling efficiency is mostly unaffected by the porosity 
and remains nearly constant, but the stress-based figure of merit 
is not diminished with the reduced compliance. It is interesting 
to compare the simulation results presented here to the poling- 
piezoelectric simulations presented by Roscow and co-workers 
[40]. For porous barium titanate simulated from zero to 40% poros- 
ity (accounting for poling and electromechanical properties), d33 
was reduced by about a half, from ~140 to ~ 80 pm/V, the elec- 
tromechanical coupling efficiency (k?) decreased from 0.45 to 0.25, 
and FoMj, was found to increase from ~1.5 to ~1.9 1071? m?/N. 
These results are a manifestation of the competing mechanisms 
discussed above. 

In comparison, for the same range of porosity in GaN, d33 in- 
creased from 3.4 to 8 pm/V (remarkably close to the experimental 
value), k? was slightly increased, and FoMį; increased 10 fold from 


~0.12 to ~1.2 107!?m?/N - close to the barium titanate. This is a 
striking result, indicating that for stress-based applications porous 
GaN can perform almost as well as barium titanate - a material 
with a piezoelectric charge coefficient (d) about 50 times higher, 
when comparing bulk forms. 


5. Conclusions 


We studied the effect of porosity on the piezoelectric proper- 
ties of a layer of conductive GaN. 40% volume porosification was 
found to deplete the structure, resulting in the onset of piezoelec- 
tricity. This was characterized at the nanoscale by PFM, where a 
piezoelectric coefficient 2-3 times larger than bulk GaN was found 
(8 compared to 2.5 - 3.5 pm/V). Macroscale energy harvesting ex- 
periments demonstrated that the porous GaN layer operates as 
a piezoelectric transducer, with maximal power densities of 100 
nW/cm/? transferred to a resistive load of ~10MQ, under 150 kPa 
load. 

We have performed detailed finite element simulations to ex- 
amine the effect of porosity on the piezoelectric coefficient and ef- 
ficiencies. The simulations suggest that the electromechanical cou- 
pling coefficient (k?) remains nearly constant with porosity - indi- 
cating it is determined by intrinsic material properties. Conversely, 
the stress-based electromechanical figure of merit, which is inde- 
pendent of effective compliance, is increased 10 fold for 40% poros- 
ity, bringing the porous GaN to a level comparable to piezoelectric 
ceramic such as barium titanate. 

Our combined modelling/experimental study is promising when 
considering porous GaN as an efficient piezoelectric material, and 
it opens up interesting avenues for integration of electromechan- 
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ical applications within dominant microelectronic processes and 
technologies. 


6. Materials and methods 
6.1. GaN growth and porosification 


n-doped GaN was grown on sapphire by metal-organic vapor 
phase epitaxy (MOVPE) and porosified in a fashion similar to pre- 
vious studies [18,19]. Briefly, a layer of n-doped GaN (~ 5:1018 
cm~?; 1.5 pm) was grown on top of a non-intentionally doped low- 
dislocation density GaN pseudo-substrate [53]. The sample was im- 
mersed in oxalic acid (0.25 M) with the nitride surface of the 
doped layer exposed and a DC potential of 8-14 V applied between 
the sample (thus determining the level of porosification) and an 
inert platinum counter electrode. The sample used in AFM and en- 
ergy harvesting measurements were etched using 14 V bias, for a 
duration such that the 1.5 um doped layer was etched to about 
a fifth of its thickness (see Supporting Information S1 and Figure 
S1). The porosification process does not affect the composition or 
crystallinity of the remaining material [18,19]. 


6.2. Atomic force microscopy (AFM) 


AFM measurements were carried out using a Bruker Dimension 
Icon, and the results were analyzed using Nanoscope Analysis 1.9 
(Bruker). For PFM and KPFM MESP-RC-V2 tips were used while for 
quantitative nanoscale mechanical mapping (QNM) the peak-force 
mode standard ScanAsyst-Air tip was used. PFM was performed by 
applying different AC voltage amplitudes at a frequency of 125 kHz 
which is away from any contact resonance. This ScanAsyst-Air tip 
is not suitable for quantitative mechanical analysis of materials as 
hard as GaN, leading to the QNM results treated as qualitative (see 
Supporting Information S2). 


6.3. Device fabrication and energy harvesting 


The device was fabricated as follows. The sample was mounted 
on a round (15 mm in diameter) metallic stub (used initially for 
AFM - Fig. 1), which was electrically connected to the conductive 
GaN layer. This served as a bottom electrode. The top electrode 
was sputtered with gold (~80 nm), made on the etched part of the 
sample using a shadow mask with a size of about 0.25 x 0.5 cm? 
(Fig. 3). 

Short-circuit and loaded mechanical energy harvesting experi- 
ments were carried out using a linear motor setup, set to tap the 
device (max. 1.5 N at 2 Hz). The generated current was measured 
using a picoammeter (Keithley 6487). The device was mounted on 
a force sensor (AEP Transducers TCA type, s-shaped), and the force 
was recorded separately from the electrical output. 


6.4. COMSOL multiphysics 


COMSOL Multiphysics v. 5.2 was used to simulate the piezo- 
electricity of the porous GaN relying on the Piezoelectricity mod- 
ule. GaN was used with an addition to a user-defined d matrix, as 
discussed above, since GaN is not in the default library of piezo- 
electric materials. The simulation was run calculating stress-charge 
quantities, and the electrical connectivity was open-circuit, hence g 
quantities were directly extracted from simulation results, yielding 
d as described above. 
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